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ABSTRACT 

In this work, we focus on the analysis of a CME-driven shock observed by SOHO/LASCO. 
We show that white-light coronagraphic images can be employed to estimate the compression 
ratio X - pdlpu all along the front of CME-driven shocks. X increases from the shock flanks 
(where X - 1.2) to the shock center (where X ^ 3.0 is maximum). From the estimated X 
values, we infer the Alfven Mach number for the general case of an oblique shock. It turns out 
that only a small region around the shock center is supercritical at earlier times, while higher 
up in the corona the whole shock becomes subcritical. This suggests that CME-driven shocks 
could be efficient particle accelerators at the initiation phases of the event, while at later times 
they progressively loose energy, also losing their capability to accelerate high energy particles. 
This result has important impUcations on the localization of particle acceleration sites and in 
the context of predictive space weather studies. 

Subject headings: Sun: coronal mass ejections (CMEs); shock waves; 

1. Introduction 

Solar energetic particle (SEP) events represent one of the most severe hazards in the near-Earth space 
environment. The particle population of SEP events at Earth critically depends on the location of the particle 
acceleration, which may occur in both the solar corona and interplanetary space. While the association 
between shocks driven by Coronal Mass Ejections (CMEs) and gradual SEP events is clear, the acceleration 
mechanisms underlying this association is less understood (Zank et al. 2000). In general, shocks take 
the flow energy upstream and convert it to thermal energy downstream and accelerated particles. Only 
supercritical shocks, i.e., shocks with Alfenic Mach number Ma (the ratio of the upstream flow speed along 
the shock normal to the upstream Alfven speed) > M* (the first critical Mach number) are thought to be 
able to reflect ions and thus produce Sjpg close to the Sun (< 3 R©). In fact, above M^, resistivity alone is 
unable to provide all the dissipation needed for the required Rankine-Hugoniot shock jump (e.g., Edmiston 
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& Kennel 1984). Critical Mach numbers range between 1 and 2 for quasi-parallel shocks and between 2.3 

and 2.7 for quasi-perpendicular shocks for the coronal plasma, depending on both Bsn and the plasma yS of 
the ambient medium through which the shock propagates. As supercriticality is surpassed, ion reflection 
is required to provide the required downstream heating (e.g., Burgess 2006). Supercritical shocks are also 
thought to be required for accelerating electrons in the energy range relevant to the production of type II 
bursts, which are the radio signature of coronal shock waves (Gopalswamy et al. 2010). Large increases 
in electron fluxes in the energy range 2-20 keV, associated with interplanetary shocks with high Mach 
number have been actually reported (e.g., Tsurutani & Lin 1985) and the fact that type II radio bursts also 
show considerable variability may be due to the associated shock changing Mach number as it propagates 
through an inhomogeneous solar wind (e.g.. Knock et al. 2003). Thus far, however, the clear identification 
of the shock-associated particle acceleration region in the corona and its evolution with time is still missing, 
although some progress has been made in the identification of candidate shock signatures both in whit- light 
images (e.g., Ontiveros & Vourlidas 2009; Liu et al. 2009) and in ultraviolet spectra (e.g., Raymond et al. 
2000; Ciaravella et al. 2005; Mancuso & Avetta 2008; Bemporad & Mancuso 2010; Mancuso 201 1). 

In this Letter, we analyze white-light observations of a fast CME-driven shock observed by the LASCO 
instrument on board SOHO. From the observed compression ratio along the shock front, we estimate M\ 
as a function of latitude and height in the corona, showing that the shock actually becomes supercritical, 
and thus more efficient in accelerating both electrons and ions to high energies and, correspondingly, in 
generating high fluxes of energetic particles, only at selected angles around the shock center. This region 
is moreover seen to disappear as the shock propagates in the coronal environment, thus implying a lower 
ability in accelerating Sjpg at progressive heights. 

2. Observations and data analysis 

Figure 1 shows SOHO/L ASCO C2 base diflFerence images of the 1999 June 11 event. This CME was 
first detected at 1 1 :26 UT in the LASCO C2 field of view below a height of 2.6 R©, travelling in the Northeast 
direction at a speed of 1569 km s~\ as provided by the online LASCO CME catalog (Gopalswamy et al. 
2009). According to GOES data, the soft-X ray emission of the associated C8.8 flare started at 11:07 UT 
and reached a maximum at 11:57 UT. The flare was emitted from active region NOAA 8585, which was 
located at a latitude of 38° above the equator and just a few degrees behind the Umb. Radio images acquired 
~ 6 minutes after the flare onset by the Nangay Radioheliograph (NRH, RadioheUograph Group 1989) show 
the formation of an arch-shaped radio emission above the limb at the latitude of the CME source region. 
The 164 MHz emission detected by NRH, which is simultaneous to the type II radio signature observed 
in ground-based metric dynamic spectra, is strongly suggestive of the plane-of-the-sky section of a CME- 
driven coronal bow-shock surface. 

The CME front was observed in only two LASCO C2 frames acquired at 11:26 and 11:50 UT. The 
derivation of the electron density from the polarized brightness (pB), which is photospheric white-light 
scattered off coronal electrons, is a classical problem (Van de Hulst 1950). By assuming spherical symmetry, 
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Fig. 1. — Base difference LASC0/C2 images showing the location of the CME-driven shock front (dashed 
lines) at 11:26 UT (left) and 11:50 UT (right). 



we fitted two-component radial power laws to the observed pre-CME pB radial profiles, from which we 
obtained a set of coronal electron density profiles all over the region of shock propagation. In order to 
determine the location of the shock front, identified as a weak white-light intensity increase located above 
the expanding CME front (Figure 1, dashed lines), we subtracted the average pre-CME intensity from each 
frame (see also Ontiveros & Vourlidas 2009 for a discussion on white-light shock identification). For each 
pixel along the front, we extracted the unpolarized white-light brightness of the shocked plasma, ^post> 
while the last image acquired before the CME arrival provided the pre-shock brightness, bp^e, at the same 
locations in the corona. In each pixel, the brightness b(() observed at a projected altitude ^ is given by a line- 
of-sight (LOS) integration of the inferred electron density profile, n^, multiplied by a geometrical function. 
The integration along the LOS is then divided into two integrals: one performed over the unshocked corona 
(with density n^) and the other over a length L across the shocked plasma with density X-ng. Here, X = pd/pu 
is the unknown shock compression ratio and pd (pu) is the downstream (upstream) plasma density. In this 
work L has been estimated as in Bemporad & Mancuso (2010), i.e. by assuming that the shock surface 
has the 3-D shape of an hemispherical shell with thickness equal to that observed on the plane of the sky, 
corrected for the shock motion during the LASC0/C2 exposure time. Resulting values are L = 0.28 R© and 
L = 0.61 Rq for the shock observed at 11:26 and 11:50 UT, respectively. Given L, by adopting the radial 
density profile derived from the analysis of the pB, the shock compression ratio X can be inferred from the 
observed tpj-^ and bpo^t. 

The so-obtained compression ratios X along the shock fronts at 11:26 and 11:50 UT, are shown in 
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Fig. 2. — Left. Compression ratios X = pd/Pu as measured at the points illustrated by the dashed lines in 
Figure 1 along the shock fronts at 11:26 UT (top) and 11:50 UT (bottom). Middle: Alfven Mach numbers 
Ma for perpendicular {dashed line) and parallel (dotted line) and for angles measured along the actual 
shock front (solid red lines) at 11:26 UT (top) and 11:50 UT (bottom); the latter curves are compared with 
the corresponding critical Mach numbers (solid blue lines). Right: Altitudes of points located at different 
latitudes along the shock front in the previous plots. 



Figure 2 (left panels). The X values vary between 1.2 and 3.0, a range that is lower than the upper limit 
compression of 4 given hy X = (y + \)l(y - 1) for a mono-atomic gas (y - 5/3). An interesting new 
result is given by the latitudinal dependence: at both times X maximizes at the center of the shock surface, 
progressively decreasing towards the flanks of the shock. Notice that values plotted in this Figure may also 
depend on unknown latitudinal changes in the L values, which have been assumed constant at all latitudes. 
A second new result is the time evolution: as the shock expands, the latitudinal dependence is preserved, 
but the X values decrease all along the shock surface. This means that the shock is losing its energy with 
propagation and that, at least for this event, it is stronger at the center of the front, as one would expect 
because this is the shock part moving faster in the corona, hence with larger Ma. The latter conclusion is 
demonstrated in the next section. 

The relationship between the compression ratio X and the Alfvenic Mach number Ma depends on 
the plasma-to-magnetic pressure ratio yS and the angle Osn between the shock normal and the upstream 
magnetic field. In the solar corona, the plasma is mainly controlled by the magnetic field (jS «; 1). In 
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the limiting case of yS ^ 0, the Alfvenic Mach number is given by Max = '^X{X + 5)/ [2(4 - X)] for a 
perpendicular shock and Ma|| - for a parallel shock (e.g., Vrsnak et al. 2002). The angles Ob,, have been 
estimated pixel by pixel along the shock surfaces from the LASCO images by assuming a radial magnetic 
field above ~ 2 R©. By taking, as a first order approximation for the more general case of an oblique 
shock, Ma/ - ^J{Mfi^j_ sin Obh)^ + {Ma\\ cos Obh)^, implying Max ^ Maz > Ma||, we finally obtained the 
values plotted in Figure 2 (middle panels). Ma/ is seen to maximizes at the shock center at both times, 
reaching a maximum value of Maz - 1.8 at ~ 2.6 Rq. Critical values for the Alfvenic Mach number M* 
for coUisionless non-relativistic shocks are provided by Edminston & Kennel (1984), recently reviewed by 
Treumann (2009). The above authors provide a theoretical M* curve for a fast-mode shock as a function of 
the angle Osn for p = 0. From the Bsn values derived along the shock front, we estimated M* at different 
latitudes. The resulting M* variations along the front are seen to be opposite to those of Ma/. That is, M* 
is seen to minimize at the shock center and maximize at the flanks. Comparison between the Maz and M* 
curves in Figure 2 shows that, at earlier times and lower altitudes, Maz > M* at the center of the shock and 
Maz < M* at the shock flanks, while, at later times and higher altitudes, Maz < M* at all latitudes. Hence, 
at earlier times the shock center is super-critical, making it a probable source for SEP acceleration. 




Fig. 3. — Cartoon showing the possible evolution of supercritical (red) and subcritical (blue) regions over 
the shock surface. 
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3. Discussion and conclusions 

In this work, we have shown for the first time that coronagraphic wliite-hght images can be employed to 
estimate the shock compression ratio X not only at a single latitude and altitude (as in previous works), but all 
over the observed shock front. The compression ratios computed for an event observed by SOHO/LASCO- 
C2 have been used to estimate the Alfvenic Mach numbers for perpendicular, parallel and oblique shocks. 
The results show that at earlier times the CME-driven shock is supercritical at the top of the shock and 
subcritical at the flanks, with the supercritical region disappearing with time. The above results suggest 
the overall evolution qualitatively shown in Figure 3, with the CME-driven shock primarily super-critical 
in its earUer phase (and hence an efficient particle accelerator). In the following minutes, as the shock 
progressively lose its energy, the supercritical region becomes narrower and narrower, with the particle 
acceleration site progressively more localized in space. At the end of the process, the whole shock becomes 
subcritical, so that no more efficient high energy particle acceleration is possible. While the present result 
for a single event supports SEP generation at the CME/shock front in the higher corona, type II emission at 
the CME flanks is still expected in the lower corona due to interaction with nearby streamers (e.g., Mancuso 
& Raymond 2004), so that a more extended study will be required in order to draw a definite conclusion. 
These results have very important impUcations on the localization of particle acceleration sites during CMEs 
and on the temporal evolution of SEP fluxes in CME-driven shocks. 

A.B. acknowledges support from ASI/INAF 1/023/09/0 contract. SOHO is a project of international 
cooperation between ESA and NASA. 
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